function l = l_i(pm,x,e,p,w,r)

l = (pm.rho./pm.theta(r).*p./w.*x.*e).^(1./(pm.theta(r)-pm.rho));